New Sensitivity to Solar WIMP Annihilation using Low-Energy Neutrinos 
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Dark matter particles captured by the Sun through scattering may annihilate and produce neutri- 
nos, which escape. Current searches are for the few high-energy neutrinos produced in the prompt 
decays of some final states. We show that interactions in the solar medium lead to a large number 
of pions for nearly all final states. Positive pions and muons decay at rest, producing low-energy 
neutrinos with known spectra, including v e through neutrino mixing. We demonstrate that Super- 
Kamiokande can thereby provide a new probe of the spin-dependent WIMP-proton cross section. 
Compared to other methods, the sensitivity is competitive and the uncertainties are complementary. 



PACS numbers: 95.35.+d, 95.85.Ry 

Introduction. — If dark matter is a thermal relic of 
the early universe, then its self-annihilation cross section 
is revealed by its present mass density. To match obser- 
vations, the required cross section, averaged over relative 
velocities, is (ctav) — (5.2 — 2.2) x 10~ 26 cm 3 /s, as a func- 
tion of increasing mass, m x TJ. This indicates a weakly 
interacting massive particle (WIMP; denoted x) 0t3|- 

The total annihilation cross section, including all fi- 
nal states, is well defined, but the partial annihilation, 
scattering, and production cross sections with any spe- 
cific standard model (SM) particles are model-dependent. 
Measurement of any of these cross sections would dra- 
matically constrain WIMP models and eliminate more 
exotic possibilities. There are limits from straightforward 
searches for astrophysical fluxes of annihilation products, 
direct nuclear scatterings in underground experiments, 
and collider events with missing energy. 

Searches for high-energy neutrinos from the Sun give 
strong limits on WIMP-nucleon scattering. WIMPs pass- 
ing through the Sun may rarely scatter with nuclei 
and become gravitationally bound. Scattering can oc- 
cur by coherent spin-independent (SI) or valence spin- 
dependent (SD) interactions. Further scatterings ther- 
malize WIMPs in the solar core, where they annihilate. 
Only neutrinos can escape and potentially be detected. 
When the capture rate, Tc, regulates the annihilation 
rate, Ta, as expected, an upper limit on the neutrino 
flux sets an upper limit on the WIMP-nucleon scattering 
cross section. 

The most interesting limits using searches for high- 
energy neutrinos from the Sun are on SD WIMP-proton 
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Though the searches are based on 



the annihilation process, these limits are independent of 
{uav), except for assumptions about the annihilation fi- 
nal states, which govern the detectability of the neutri- 
nos. High-energy neutrinos come from the few annihi- 
lation products that decay promptly, before losing en- 
ergy in the solar medium, giving continuum spectra up to 
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(direct annihilation to neutrino pairs, helicity- 



suppressed in many models, produces a line at E v = m x ). 
Strong limits on the SD WIMP-proton scattering cross 
section have been derived for large m x and certain final 
states, such as W + W~ ,t + t~ , and bb. 

Generalizing to other final states and a broader range 
of masses is challenging. The number of high-energy neu- 
trinos per annihilation is small (N u ~ 1 for favorable final 
states) and their spectrum depends on the unknown final 
state. Neutrino signal detection and atmospheric neu- 
trino background rejection become easier for larger m x , 
though neutrinos with E v > 100 GeV are significantly 
attenuated by interactions in the Sun. 

What about the more common but seemingly less fa- 
vorable final states? As is well known, most final states 
ultimately produce pions and muons, which quickly lose 
energy and decay at rest, producing only MeV neutrinos, 
which have long been considered undetectable (e.g., see 
Ref. [B|, Q and many subsequent papers) . 

We propose a new probe of SD WIMP-proton scat- 
tering in standard solar WIMP capture scenarios. First, 
we show that the pion yield from WIMP annihilation in 
the Sun - produced directly in hadronic decays and fur- 
ther through inelastic interactions in the dense medium 
- is large, relatively model-independent, and increases 
as N-x oc m x . Second, we show that the subsequent 
low-energy neutrinos are much more detectable than 
previously thought, due to their high yield (AT„ 3> 1) 
and known spectra, the low atmospheric neutrino back- 
grounds at low energy, and advances in detectors. The 
advantages of this method are sensitivity to low WIMP 
masses and near-inscnsitivity to the choice of final state. 

In the following, we review solar WIMP capture and 
annihilation, calculate the pion and neutrino yields per 
annihilation, estimate the signal and background rates 
in Super-Kamiokande (Super-K) [7j, and derive new con- 
straints on the SD WIMP-proton cross section. We con- 
clude by discussing likely improvements and the impor- 
tance of complementary methods to allow for possible 
surprises in the astrophysics or physics of dark matter. 
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FIG. 1: Solar WIMP capture rate as a function of WIMP 
mass for two values of the SD WIMP-proton scattering cross 
section. Evaporation (see text) is not included. 

WIMP Capture and Annihilation. — When the 
probability of a WIMP scattering with a nucleus while 
passing through the Sun is small, the capture rate scales 
linearly with (we consider just the dominant process 
of SD WIMP-proton scattering). It also scales linearly 
with the number density, p x /m x , where p x is the local 
mass density of WIMPs. For large masses, the capture 
rate falls more rapidly than l/m x , due to kinematic sup- 
pression of the energy loss [1, [j| . 

Figure Q] shows r c computed with DarkSUSY 0. We 
use the defaults of p x = 0.3 GeV/cm 3 , a Maxwellian ve- 
locity distribution for WIMPs with a 3-D velocity dis- 
persion of 270 km/s, and a circular velocity of the Sun of 
220 km/s. For the parameters considered here, the anni- 
hilation rate easily reaches equilibrium with the capture 
rate, so Ta = Fc/2. We use DarkSUSY only to calcu- 
late the capture rate. We treat the annihilation processes 
through a dedicated separate simulation, as described be- 
low. Scattering interactions in the hot solar core can un- 
bind WIMPs ( "evaporation" ) ll| , and this process is not 
included in DarkSUSY. Detailed calculations show that 
it is important for low WIMP masses but negligible for 
larger masses. To be conservative, we present our final 
results for only m x > 4 GeV, though less stringent limits 
could be calculated for lower masses. 

Pion Yield from Annihilation. — We calculate the 
number of pions produced per annihilation as a function 
of WIMP mass. This depends on (i) the fraction of the 
annihilation energy in final states with eventual hadronic 
content, (ii) the hadronization of quarks and gluons, and 
(iii) all decay, interaction, and energy-loss processes in 
the solar medium. The first is easily parameterized, the 
second is the same as in vacuum, but the third requires 
a detailed simulation. 

For hadronic final states, we adopt results measured 
at e + e _ colliders 12j: the initial pion multiplicity from 
hadronization is iV* ~ 3+4.51og^ (y / s/GeV), pions dom- 
inate over other hadrons by an order of magnitude, and 



FIG. 2: Low-energy neutrino yield (summed over flavors) 
per WIMP annihilation (y/s = 2m x ) in the solar core, ob- 
tained by simulating pion-induced hadronic showers in the 
solar medium. Reference lines are shown for the cases where 
all or only 10% of the annihilation energy goes into producing 
pions. The high-energy neutrino yield for r + r~, relevant for 
current searches, is also shown. 

there are comparable populations of tt + , 7r~, tt°. In vac- 
uum, hadrons simply decay, producing a modest number 
of additional pions. In matter, by contrast, the number 
of pions can be greatly amplified. 

For high-energy charged pions, the hadronic interac- 
tion length in the Sun is short compared to the decay 
and continuous energy-loss lengths, so the number of pi- 
ons increases in each generation of the hadronic shower 
until loss processes dominate at low pion energies. Neu- 
tral pions only decay, due to their short lifetime, divert- 
ing 1/3 of the hadronic energy in each generation into an 
electromagnetic shower. Charged pions are eventually 
brought to rest. Negative pions are Coulomb-captured 
into pionic atoms with nuclei beyond hydrogen and then 
absorbed, and so do not produce neutrinos [1 31 ] . Positive 
pions decay, producing 3 neutrinos with energies up to 
E v = 52.8 MeV, through tt + -> [i + v^ -> e+VeV^. 

We used GEANT4 [3| to model hadronic showers in 
the solar core, defined by a homogeneous volume with 
density 160 g/cm 3 and solar elemental abundances fl5j |. 
Each annihilation releases \fs = 2m x of energy, which 
we assume goes to a hadronic final state with the ini- 
tial pion multiplicity and equal numbers of ir + , ir~ , 
7T°. We inject these into the solar core and calculate the 
resulting final pion yield N v . We find that our results 
are not sensitive to reasonable variations in the density, 
composition, or ionization state. At these energies, pion 
energy losses on photon targets are negligible [1 61 ) - 

Figure [2] shows our results for the low-energy neutrino 
yield, N v — 3N V +, per WIMP annihilation. The near- 
linear scaling of N v with annihilation energy indicates the 
primary importance of the total hadronic energy, with 
the mild deviations reflecting losses due to ir° produc- 
tion. Comparable results are indeed obtained when the 
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initial number of pions is varied at fixed injection energy. 
Because of pion dominance in hadronic showers, similar 
results are expected for arbitrary hadronic annihilation 
products, and this was confirmed in some representative 
cases. Our calculated results are comparable to simply 
assuming that ~ 10% of the annihilation energy goes 
into producing pions, i.e., N„ ~ 0.1 (2m x /m^). This 
is similar to the fraction of energy in hadronic showers 
in ultra-high-energy cosmic ray interactions in Earth's 
atmosphere [T3|; the details differ, but the common el- 
ement is that the dominant charged-pion loss process is 
hadronic scattering. 

For low-energy neutrinos, the rate and spectrum are 
very similar for nearly all final states, unlike the case for 
high-energy neutrinos. Quarks, gluons, weak bosons, and 
tau leptons all decay dominantly into hadrons, which pro- 
duce pions. Muons also produce low-energy neutrinos. 
High-energy neutrinos are covered in current searches. 
For gamma rays and electrons, there are very few neutri- 
nos produced (however, those final states are well probed 
by searches for annihilation fluxes from Galactic and ex- 
tragalactic dark matter structures). 

Neutrino Signals and Backgrounds. — The neu- 
trino flux is 0„ = N U T a/{^D 2 ), where D is the Earth- 
Sun distance. The neutrino spectra are well known: 
7T + — > gives a line spectrum for i/^, and /x + — ¥ 

e + u e u^ gives a softer continuum spectrum for v e and a 
harder one for fl2| . The Sun and Earth are transpar- 
ent to these neutrinos. We consider detection in Super- 
K, using the interaction v e p — > e + n. Matter-enhanced 
mixing in the Sun gives a flavor-change probability of 
P{y^ — » D e ) ~ 1/6 for either neutrino mass hierarchy [18]. 

The detectability considerations are similar to those 
for the diffuse supernova neutrino background, though 
that peaks at lower energies and monopole- 

catalyzed proton decays in the Sun, though that has 
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lower pion injection energies [21 
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For this neutrino 
cross section rises as a v ~ 
10~ 43 [E v j MeV) 2 cm 2 and the positron is detectable; 
it has nearly the full neutrino energy, but only a weak 
forward anisotropy in this energy range. Super-K has 
1.5 x 10 33 free (hydrogen) protons in the 22.5 kton of 
water in the fiducial volume. 

Figure [3] shows the observable v e signal spectrum, ac- 
counting for neutrino mixing, the full neutrino cross sec- 
tion and kinematics, and detector e nerg y resolution Q; 
the detection efficiency is near unity [20|. For simplicity, 
we consider just the 4.1 livetime years of the Super-K-I 
phase (1996-2001). The measured background spectrum, 
which has a small absolute rate, is also shown; its compo- 
nents are displayed in Ref. [2(J. The largest component 
is electrons and positrons from the at-rest decays of sub- 
Cerenkov ("invisible") muons produced by atmospheric 
neutrinos. The spectrum of this background is the same 
as the signal (which mixes to v e ) flux spectrum. How- 
ever, the detection spectrum of the signal is weighted 
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FIG. 3: Detectable signal in Super-K induced by low-energy 
v e from solar WIMP annihilation, along with the measured 
background (both 4.1 livetime years). The signal shape is 
independent of WIMP properties, and its normalization scales 
with er^p (here chosen to be at the edge of exclusion). 



with the neutrino interaction cross section. 

Constraints on WIMP Properties. — Figure 0] 
shows the Super-K sensitivity to er^ with present data. 
We estimate this by first calculating the expected num- 
bers of signal and background events, summed over the 
energy range 25-50 MeV. These numbers are 15 and 135 
for the case shown in Fig. [3] For each WIMP mass, we 
derive the 90% C.L. upper limit on that would be 
set if no signal were seen and statistical fluctuations in 
the background were taken into account. Our approach 
is conservative, and improvements are discussed below. 

Our estimate is compared with published upper limits 
from direct nuclear scattering experiments (see also new 
results in Ref. (24|) and indirect searches for high-energy 
neutrinos from the Sun. The sensitivities of both are 
expected to improve greatly with better detectors and 
analyses. However, it is very difficult for each to ex- 
tend to lower m x , as the signals decrease while the back- 
grounds increase, and there are kinematic thresholds. In 
contrast, the sensitivity of our proposed low-energy neu- 
trino signal improves with decreasing m x . There are also 
strong limits on the SD WIMP-nucleon scattering cross 
section that can be deduced from limits on mono-jet and 
mono-photon signals at hadron colliders; these require 
that the masses of mediators coupling WIMPs to the SM 
be large [25ll26j. 

Conclusions. — We propose a new probe of the SD 
WIMP-nucleon scattering cross section, using low-energy 
neutrinos produced through pion multiplication and de- 
cay following WIMP annihilations in the Sun. We es- 
timate the prospects for Super-K, finding sensitivity to 
cr^p in a range competitive to that of other, rather differ- 
ent, experiments. Importantly, our results are nearly in- 
sensitive to the annihilation final states and the details of 
the astrophysical inputs (in standard scenarios [27j); non- 
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FIG. 4: New sensitivity to using low-energy v e in Super- 
K. Possible improvements are discussed in the text. For m x < 
4 GeV, evaporation of solar WIMPs will degrade the si gnal 



Published upper limits (90% C.L.) from direct scattering 
[30l ] and high-energy neutrino searches are shown [3~fl. 
along with the possible DAMA/LIBRA signal region [33| 
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tions through a new low-mass force carrier that decays 
only into low-mass SM particles 41] could have vanish- 
ing high-energy but strong low-energy neutrino signals. 
These are just a few examples. Most generally, a com- 
bination of experiments will be required to reconstruct 
or constrain the complete couplings of the SM to dark 
matter, whether WIMPs or something more exotic. 

Note added: As this work was being completed, we 
learned of an independent study with similar results by 
Bernal, Martin- Albo, and Palomares-Ruiz [42[ that will 
be submitted to arXiv simultaneously. 
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standard scenarios can be different and are mentioned 
below. In addition, the sensitivity easily extends to the 
region of low masses, which is currently of great interest 
and is challenging to probe with other methods. A ded- 
icated study by the Super-K collaboration, using more 
data, neutrino-electron scattering signals for all flavors, 
and full energy spectra and angular distributions, should 
give immediate improvements on our estimates. 

In the future, Super-K may be enhanced with dissolved 
gadolinium to allow neutron detection and thereby bet- 
ter separation of signals and backgrounds [35| . If a com- 
bination of techniques removed the backgrounds, then 
the sensitivity to the neutrino signal and hence also cr^p 
could improve by ~ 15, beyond which even Super-K is 
too small to expect any signal events. Proposed large liq- 
uid scintillator [36[ or liquid argon [37} detectors would 
also have interesting sensitivity. Hyper-Kamiokande [38[ , 
which is intended to be 25 times larger than Super-K 
and which also may have gadolinium, could potentially 
improve on present estimates by up to ~ 15 x 25 ~ 375. 

Direct, indirect, and collider probes of the SD WIMP- 
proton scattering cross section rely on different assump- 
tions and hence are complementary. Multiple methods 
with comparable sensitivity are needed to test results 
from one search against the others. These tests could 
provide deep insights into the astrophysical distributions 
and particle properties of WIMPs. 

If common assumptions about dark matter are incor- 
rect, then the relative power of different methods could 
change dramatically. A dark matter disk with low- 
velocity WIMPs could enhance both types of neutrino 
signals [3^|. A time- varying WIMP flux from dark mat- 
ter substructures could alter the relationship between 
neutrino and scattering signals [Io[. WIMP annihila- 
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